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Abstract 

MAGIC is a system of two Imaging Atmospheric Cherenkov Telescopes located in the Canary island of La Palma. Since autumn 
2009 both telescopes have been working together in stereoscopic mode, providing a significant improvement with respect to the 
previous single-telescope observations. We use observations of the Crab Nebula taken at low zenith angles to assess the performance 
of the MAGIC stereo system. The trigger threshold of the MAGIC telescopes is 50 - 60 GeV. Advanced stereo analysis techniques 
allow MAGIC to achieve a sensitivity as good as (0.76 + 0.03) % of the Crab Nebula flux in 50 h of observations above 290 GeV. 
The angular resolution at those energies is better than ~ 0.07°. We also perform a detailed study of possible systematic effects 
which may influence the analysis of the data taken with the MAGIC telescopes. 
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1. Introduction 

MAGIC (Major Atmospheric Gamma Imaging Cherenkov 
telescopes) is a system of two 17 m diameter Imaging Atmo- 
spheric Cherenkov Telescopes (IACT). They are located at a 
height of 2200 m a.s.l. on the Roque de los Muchachos on the 
Canary Island of La Palma (28°N, 18°W). They are used for 
observations of particle showers produced by very high energy 
(VHE, > 30 GeV) gamma-rays from galactic and extragalactic 
objects. 

While the first MAGIC telescope has been operating since 
2004, the construction and the commissioning of the second 
telescope finished in autumn 2009. Both telescopes are nor- 
mally operated together in the so-called stereoscopic mode. 
In this mode, only events which trigger both telescopes are 
recorded and analyzed. Since the same event is seen by two 
telescopes, with a slightly different parallax angle, it allows for 
the full three-dimensional reconstruction of air showers. 

The Crab Nebula is a nearby (~ 2 kpc away) pulsar wind 
nebula and the firs t source detected in VHE gamma-rays 
dWeekes et allll989h . Moreover, the Crab Nebula is the bright- 
est steady VHE gamma-ray source, therefore it has become 
the so-called "standard candle" in VHE gamma-ray astronomy. 
Recent obse r vation of flares in the G eV range (ATel #2855, 
Tavani et all 1201 lc lAbdo et all 1201 II) raised questions about 



system is referred to as the active mirror control (AMC). The 
light mechanical structure of both telescopes allows for a rapid 
repositioning time, necessary for observations of short phenom- 
ena such as gamma-ray bursts. The tel escopes can perfor m a 



the stability of the TeV flux. In fact a hint of increased flux was 
claimed by ARGO-YBJ (ATel #2921), however it was not con- 
firmed by the simultaneous observations of MAGIC and VERI- 
TAS (ATels #2967, #2968). In this paper we use the Crab Neb- 
ula data to evaluate the performance of the MAGIC telescopes. 
The scientific implications of the Crab Nebula observations and 
the stability of its flux will be addressed in a separate paper in 
preparation. 

In Section [2] we briefly describe the MAGIC telescopes. In 
Section [3] we describe both, the Crab Nebula data sample, and 
the various Monte Carlo samples used in the analysis. In Sec- 
tion|4]we present the techniques and methods used for the anal- 
ysis of the MAGIC stereo data. In Section [5] we evaluate the 
performance of the MAGIC stereo system. In Section [6] we 
analyze and quantify various types of systematic errors. We 
conclude the paper with final remarks in Section|7] 

2. The MAGIC telescopes 

The second MAG IC telescope is a close copy of MAGIC I 
( Albert et al. , 2008bl) with a few important differences. Both 
telescopes have a 17 m diameter mirror dish and f/D close to 
1, but contrary to MAGIC I, the outer mirrors of MAGIC II are 
made of glass rather than aluminium dDoro et all 120081) . Due 
to the parabolic shape of the reflector, the time spread of syn- 
chronous light signals is very low (< 1 ns). To correct for de- 
formations of the dish and sagging of the camera during the ob- 
servations, each mirror panel is equipped with 2 motors. This 



half-turn of 180° in azimuth in just 20 s dBretz et alH2009t) 



In order to be able to register short light pulses, each MAGIC 
telescope is equipped with a 3.5° diameter camera with photo- 
multipliers (PMTs) as individual pixels. In the MAGIC I cam- 
era two types of pixels are used: the inner 397 PMTs have a 
diameter of 0.1°, while the 180 outer ones are larger (0.2° di- 
ameter). The most central pixel of the MAGIC I camera has 
been modified to be able to register also constant values of the 
flowing current. This pixel is normally used only for the obser- 
vation of pulsars dLucarelli et aUl2008l) . In contrast, the camera 
of MAGIC II is homogeneous and consists of 1039 hexa gonal 



pixels with a diameter of 0.1° dBorla Tridon et all 120091) . For 
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the calibration of the PMT signals each telescope is equipped 
with a calibration box mounted in the middle of the reflector 
dish. The calibration box provides short calibration pulses of 
constant intensity uniformly illuminating the camera. The cal- 
ibration b ox of MAGIC I is base d on a set of LEDs in differ- 
ent colors dSchweizer et al. I l2002l) . while the calibration box of 
MAGIC II is composed of a frequency tripled Neodym YAG 
microchip laser and a set of attenuation filters. 

Electrical pulses at the output of the PMTs are converted 
into optical signals with the use of vertical cavity laser diodes 
(VCSELs) and transmitted over optical fibers into the count- 
ing house. There, the signal from a given PMT is split in the 
so-called receiver boards into the trigger branch and the read- 
out branch. The readout of the MAGIC I telescope is based on 
multiplexed FADCs (Flash Analog to Digital converters). The 
signals from groups of 16 pixels are artificially delayed by dif- 
ferent time offsets and are plugged into a single FADC channel, 
sampling at 2GSamples/s dGoebel et all 120071) . Although the 
same sampling frequency is also used in MAGIC II, its readout 
is based on the Domino Ring Sampler 2 (DRS2) chip, allowing 
for a lower cost and a more compact system. The readout of the 
second telescope introduces a dead time of about 10%, whereas 
the dead time of the MAGIC I readout is negligible. Ultra- 
fast PMTs and readout electronics (together with the parabolic 
shape of the reflector) are the key elements in allowing us to use 
the time information in the reconstruction of the showers. 

For each pixel in the trigger region the Individual Pixel Rate 
Control (IPRC) software controls the discriminator threshold in 
real time (the so-called level-zero trigger). It keeps the acci- 
dental event rate coming from the night sky background (NSB) 
and the electronic noise at a constant level. Each telescope sep- 
arately has a level-one trigger with the next neighbour (NN) 
topology. xNN trigger topology means that only events with 
a compact group of at least x pixels surviving the level-zero 
trigger pass the level-one trigger criterion. Unlike in mono ob- 
servations, where 4NN topology is used, stereo data require a 
3NN condition. The trigger area for MAGIC II is larger (1.25° 
radius) than for MAGIC I (0.95° radius). For special observa- 
tions aiming for the lowest possible threshold, as an alterna- 
tive to the l evel-o ne trigger, the so-called SUM-trigger is used 
dAliu et alll2008l) . Both telescopes are working in stereo mode, 
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i.e. only events triggering both telescopes are recorded. The 
level-one trigger signals of each telescope are propagated to the 
common stereo trigg er board, where the coincidence is formed 



(Pao letti et al.L 120081) . Since the telescopes are separated by 
85m the shower front will reach them at different times. Thus 
the level-one trigger signals are delayed by a value which is de- 
pendent on the pointing direction of the telescopes. With this 
procedure, the width of the coincidence gate is 100 ns. 

3. Data sample 

In order to evaluate the performance of the MAGIC Stereo 
system, we use a sample of Crab Nebula data taken in Novem- 
ber 2009 and January 2010. The data were taken at low zenith 
angles (< 30°). After a data selection based on the rate of 
background events, 9h of good quality data were obtained. 
Observations were performed in the so-called wobble mode 
dFomin et all 1 19941) . i.e. with the source position offset by 0.4° 
from the camera center in a given direction. This method al- 
lows us to simultaneously obtain the signal and background es- 
timations at identical offsets from the viewing direction of the 
telescopes. Every 20 minutes the direction of the wobble off- 
set is inverted in order to decrease the systematic uncertainties 
induced by possible exposure inhomogeneities. In addition, we 
analyze a few dedicated samples of Crab Nebula data taken at 
various offsets from the camera center. We use those samples to 
evaluate the performance of the MAGIC telescopes for sources 
at a non-standard distance from the camera center. 

For the analysis we use a sample of low zenith Monte Carlo 
(hereafter MC) gamma-rays generated with energies between 
30 Ge V and 30 Te V with the generation radius £> max = 350 m. 
In order to perform some basic background studies we also use 
diffuse MC samples of protons (30 GeV - 30 Te V, simulated on 
a square 1.2 km xl.2 km, viewcone semi-angle 6 = 5°), helium 
nuclei (70 GeV - 20 TeV, £ max = 1200 m, 6 = 6° ) and electrons 
(70 GeV - 7 TeV, Z? max = 650 m, 6 = 4.5°). 



4. Analysis 

The analysis of the data has been performed using the stan- 
dard software package called MARS (acron ym for MAGIC 



Analy sis and Reconstruction Software, see iMoralejo et al 



Analy; 
20091) . 



4.1. Image processing 

In the first steps of the analysis (cali bration, image clean - 



ing, and paramet rization of images, see lAlbert et al.l (12008bl) 



Aliu et al.1 j2009h ) of the MAGIC stereo data each telescope 
is treated separately. The procedure is similar for both tele- 
scopes with some small differences. The signal extraction in 
each channel of MAGIC I uses a cubic spline algorithm. A 
number of FADC counts (after substraction of the pedestal) in 
time slices is interpolated by a t hird degree polynomial and then 
integrated d Albert et al.Ll2008ah . In MAGIC II, the raw da ta has 



to be linearized before processing dTescaro et al. . |2009h . Af- 



terwards the signal is extracted using a sliding window with a 



width of 4 ns from the total readout window of 40 ns. The value 
of the baseline in the MAGIC II channels is estimated from the 
first 8 ns of the same readout window on an event-by-event ba- 
sis. The conversion from integrated FADC counts to photoelec- 
trons (phe) is done in both telesc opes using the F-Fac t or (ex - 
cess noise factor) method (see e.g. Mirzovan & Lorena. 1997 ). 
Due to different readout system in both telescopes, the typi- 
cal conversion factors have much different value for MAGIC II 
(~ 0.01 phe/cnts) than for MAGIC I (~ 0.002 phe/cnts for inner 
pixels). 

The time response of the DRS2 chip is not homogeneous. 
The channels of the Domino chip are subject to a time delay 
(up to a few ns) whose length depends on the absolute position 
of the readout window within the Domino ring. Using the laser 
pulses taken during calibration runs we can calibrate this effect. 
This allows us to obtain a time resolution for large, simultane- 
ous calibration pulses as good as 0.33 ns, which is of he order 
of the sampling frequency. 

A typical air shower event produces a measurable signal in 
between a few and a few tens of pixels of the camera. The 
signal obtained in the rest of the pixels is induced by the NSB 
and the electronic noise. The procedure used to select event- 
wise the pixel s which come fro m a shower is called "time im- 
age cleaning" (lAhu et all 120091) . The time image cleaning used 
in MAGIC is a two stage process that comprises searching for 
core and boundary pixels. For MAGIC I the core of the im- 
age is composed of pixels with signals above 6 phe. Pixels with 
an arrival time different by more than 4.5 ns than the mean ar- 
rival time are excluded from the core. Single, isolated pixels 
are also removed from the image in order to avoid signals from 
PMT afterpulses. An additional layer of boundary pixels with 
a signal above 3 phe and and an arrival time within 1.5 ns from 
the neighbouring core pixel is added to the image in the second 
cleaning step. For MAGIC II, due to a higher light collection 
efficiency, a larger amount of pixels, and a somewhat higher 
noise, the thresholds are increased to 9 phe and 4.5 phe for core 
and boundary pixels respectively. The time constraints are the 
same for both telescop es. Images ar e parametrized using the 
classical Hillas ellipses (IHillasL 1 1 9 8 51) . The total number of phe 
in the whole image, the S ize parameter, is strongly correlated 
with the energy of a gamma-ray event (see Fig.Q]). 

Another cleaning algorithm ; the so-calle d SUM image clean- 
ing is implemented in MARS (Rissi, l2009h . In the case of SUM 
image cleaning groups of 2, 3 or 4 neighboring pixels with sig- 
nals coming within a short time window are searched for. Each 
group is accepted if the sum of signals is above a given thresh- 
old (different for each multiplicity). Individual pixel signals 
are truncated before computing the sum to limit the effect of 
PMT afterpulses. In this paper we use the standard time image 
cleaning. The SUM image cleaning will be discussed in another 
paper (in preparation). 

4.2. Stereo parameters 

After the previous steps in the analysis chain, the im- 
ages from both telescopes are combined to obtain the three- 
dimensional stereo parameters: the Impact parameters, defined 
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Figure 1: Correlation between the mean Size of a gamma-ray event and its 
energy obtained with MC simulations. 



as the orthogonal distances of the shower axes to the tele- 
scopes, and the height of the shower maximum (hereafter the 
MaxHeight). The stereo parameters improve the energy esti- 
mation and the background rejection (lFegaiull997l) . 

The distributions of MaxHeight for different bins of image 
Size (total number of phe in the image) are shown in Fig. [2] 
The distribution of MaxHeight for gamma-rays is a Gaus- 
sian. Higher energy showers penetrate deeper into the atmo- 
sphere. Therefore, the maximum of the MaxHeight distribution 
is shifted to lower values for larger S ize. On the other hand, the 
shape of the MaxHeight distribution for background events is 
more complicated. For large images (S ize > 300 phe), hadronic 
events have a single-peaked distribution of MaxHeight, which 
is a bit broader than for gamma-ray showers. In contrast, for 
events with S ize < 300 phe a second peak appears at the height 
of ~ 2-3 km above the telescopes. This additional peak 
is pro duced by single muon events (see also Maier & Knappi 
20071) . Therefore, MaxHeight is a powerful parameter, when 
used for gamma/hadron separation at low energies. It is espe- 
cially important taking into account that small images usually 
have poorly determined Hillas parameters. 



4.3. Gamma/liadron separation 

For the gamma/hadron separation in the stereo observations 
we use the same method, the so-called Random F orest (RF), as 



has be en used for the single-telescope analysis dAlbert et al 
2008a) . The RF combines individual parameters from each 



telescope with the stereo parameters to produce a single num- 
ber, the so-called Hadronness, which characterizes the like- 
lihood of a given event to be of hadronic origin. In total, 
1 1 parameters are used for the gamma/hadron separation. In 
addition to the standard Hillas parameters of each telescope 
(log(5/zei), log(5/ze2), Widths, Vfidih-i, Length\, Lengthy) and 
the stereo parameters (Impacti, Impacts, MaxHeight) we also 
use the gradients of the arrival times of the signals in the pix- 
els projected to the main axis of the image {T imeGradient\, 
TimeGradienti). Those parameters do not assume any a priori 



known source position. Therefore, the obtained Hadronness is 
not biased by the location of the source and can also be used to 
produce a skymap of the observed region. 

The survival probability for gamma-rays and background 
events is shown in Fig.[3]for different cuts in the Hadronness 
parameter. At medium and higher energies, the Hadronness al- 
lows us to achieve a reduction of the background by a factor of 
~ 100 with only a mild loss of gamma-ray events. The bump 
and dip structure clearly visible in the Hadronness of gamma- 
rays is connected with the internal transition of the separation 
mostly based on the MaxHeight parameter (at lower energies) 
to the one based on the Hillas parameters. The separation power 
drops at very high energies, due to a high fraction of truncated 
events and small statistics in the training sample. The collection 
area of the MAGIC telescopes using cuts optimized as a func- 
tion of energy, so as to provide the best differential sensitivity 
(see Section l578l is presented in Fig. [4] Because of the threshold 
effect, the collection area after all analysis cuts is rather steep 
function of energy below 100 GeV and reaches a plateau at TeV 
energies. 

4.4. Energy reconstruction 

The energy of an event is estimated by averaging individual 
energy estimators for both telescopes. These estimators are de- 
rived from a look-up table based on a simple model for the light 
distribution of a shower based on the Size, Impact parameters 
of one telescope plus the MaxHeight parameter and the zenith 
angle of the observation. Assuming that most of the light pro- 
duced by a gamma ray in the atmosphere is contained in a light 
pool of radius r c , the mean photon density in the light pool from 
a single charged particle of the gamma-ray shower can be calcu- 
lated from the total power of emitted light by such a particle at a 
given height in the atmosphere. Given the energy, MaxHeight 
and zenith angle of a gamma-ray event, the amount of light pro- 
duced, the light pool radius (r c ) and the mean photon density in 
the pool from a single particle (p c ) are computed using a sim- 
ple atmospheric model. We use look-up tables filled with MC 
information of (E true x p c /S ize) as a function of Impact /r c and 
S ize for each telescope, to obtain the energy of an event mea- 
sured by each telescope. This parametrization takes advantage 
of the fact that E true is roughly proportional to the number of 
secondary particles in the shower maximum, which scales as 
Size/p c and that the zenith angle dependence is automatically 
included. Using the MC tables generated in the previous step, 
the energy is calculated for each telescope from the values of 
Impact, r c , p c and Size computed in each event. The final en- 
ergy estimation, E rec , is the average value obtained from both 
telescopes (weighted with the inverse of the error). 

4.5. Event arrival direction reconstruction 

Classically in the stereo systems of Cherenkov telescopes the 
arrival direction of the primary particle is estimated as the cross- 
ing point of the main axes of the Hillas ellipses in t he ind ividual 
cameras ( Aharonia n et al 1 Il997l iHofmann et all fl999h . This 
method uses only geometrical properties of the images, neglect- 
ing the timing information. 
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Figure 2: Distribution of MaxHeight in different ranges of the mean Size parameter for MC gamma-rays (dotted line), MC protons (full circles), and real data 
(solid line). 



As was shown for the case of the MAGIC I stand-alone tele- 



scope (lAleksic et all 1201 Oal) . the inclusion of timing informa- 
tion can considerably improve the angular resolution, and as a 
result also the sensitivity. In that analysis the so-called DISP 
parameter (the angular distance between center of gravity of 
the image and the estimated source position) was estimated 
using multidimensional decision trees (the so-called DISP RF 
method). Then for each image, the source position was recon- 
structed to be DISP distant from the centre of gravity along the 
main axis of the image. 

For stereo observations we apply the modified DISP RF 
method (hereafter Stereo DISP RF). First, in addition to the 
crossing point of the main axes, we calculate also the DISP 
value for images in each telescope separately. If the DISP posi- 
tions from both telescopes agree with each other within a given 
radius we compute a weighted average of them (see Fig. |5). 
However it may happen, that the direction of DISP is misre- 
constructed and the source position is located on the other side 
of the image (wrong head-tail discrimination). Therefore, if 
the reconstructed DISP positions from both telescopes are far 
apart, we check if one of the possible DISP position is close to 
the crossing point of the main axes. In the opposite case (i.e. 
none of the DISP positions agrees with the crossing point) the 
event is rejected. In this way the angular resolution is always 
obtained using the most efficient method. 

Recently, a simplified version of the above procedure, choos- 
ing the closest match among the four possible combinations of 
the two DISP positions in each telescope, was found to deliver 
a similar performance and is now used in MAGIC as the default 
algorithm. That method does not require the crossing point in- 
formation. 

In addition, since the DISP method is trained on gamma- 
rays, often both reconstructed DISP positions do not agree for 
hadrons, and the event is rejected. Therefore, an additional 
gamma/hadron separation factor is achieved at this step, which 
depending on the strength of the hadronness cut can remove be- 
tween 10% and 50% of the background events. The so-called Q 1 
plots (the distribution of the squared distance between the true 
and the reconstructed source position) in two energy ranges are 



shown in Fig. [6] Observations with the MAGIC telescopes yield 
a highly significant signal below 100 GeV. Moreover, a strong 
source such as the Crab Nebula allows for a nearly background 
free analysis above 300 GeV. 



5. Results 

5.1. Trigger threshold 

In order to achieve the lowest energy threshold, the level-zero 
trigger amplitudes for the individual pixels are chosen as low 
as possible while keeping a negligible accidental stereo trig- 
ger rate. As a result, both telescopes have a level-one trigger 
rate dominated by NSB triggers and PMT afterpulses. Typi- 
cal level-one trigger rates are around 10 kHz for MAGIC I and 
around 15 kHz for MAGIC II, due to the larger trigger area of 
MAGIC II. The resulting accidental stereo trigger rate is just 
~ 10 Hz, well below the stereo event rate from cosmic rays 
which typically is within the range 150-200 Hz. Accidental 
events are removed in the image cleaning procedure. 

The energy threshold of the stereo system is estimated from 
Monte Carlo simulations where the individual pixel thresholds 
match those used during data taking. Nominal values of the in- 
dividual pixel thresholds should be around 4 phe per 0.1° pixel, 
but the values are automatically adjusted during data taking to 
keep the accidental trigger rate at a low level. To avoid a bias 
in the conversion from discriminator threshold (amplitude mea- 
sured in mV) to photoelectrons (which is an integrated charge), 
that will translate in further systematic errors, the values of the 
discriminator thresholds are estimated from the data. This is 
done by studying the pixel amplitude of the smallest showers 
that could trigger the system. The charge distribution of the 
pixel with the lowest charge in the 3NN pixel combination with 
largest total charge (the charge of the 3 pixels added) is used to 
obtain the value of the discriminator threshold measured in phe. 
For single telescope triggers, the peak of the distribution should 
be very close to the actual number of photoelectrons needed to 
trigger a pixel. For the stereo system it is biased to higher val- 
ues and a comparison with MC simulations for different pixel 
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Figure 3: Fraction of surviving background (data, top panel) and gamma- 
ray (MC, bottom panel) events. Different line styles are for different cuts: 
Hadronness < 0.1 (solid), < 0.2 (dashed), < 0.4 (dotted), < 0.6 (dot-dashed), 
and < 0.9 (dot-dot-dashed). 



threshold settings is needed in order to find the best match with 
the data distribution. 

The distribution of the minimum single pixel charge in the 
3NN combinations computed from data taken from a galac- 
tic source observed at a low zenith angle are shown in Fig. [7] 
for MAGIC I and MAGIC EL The MC distributions also shown 
in the same figure are those that match best the measured dis- 
tribution on the data. The value of the pixel threshold is ob- 
tained from the values used in the MC simulation shown in the 
same figure. We obtain that the average threshold is 4.3 phe 
in MAGIC I and 5.0 phe in MAGIC EL From the MC simula- 
tions we also obtain that the energy threshold, defined as the 
maximum of the distribution of triggered gamma-ray events for 
sources with a differential spectral index -2.6 is 50 GeV. If we 
use instead the measured, curved Crab Nebula spectrum, the 
maximum of the distribution of triggered gamma-ray events is 
60 GeV. After applying the image cleaning and cuts used in this 
analysis, the energy threshold increases to 75 GeV. 




4 4.5 

log(E/GeV) 



Figure 4: Gamma-ray collection area of the MAGIC telescopes (MC simula- 
tions) at the trigger level (dashed line) and after all analysis cuts (image cleaning 
of both images, Size > 50 in both telescopes, gamma/hadron and 6 2 cuts). 




Figure 5: Principle of the Stereo DISP RF method. The crossing point of 
the main axes of the images is shown as an empty square, and the DISP RF 
reconstructed position from each telescope is an empty circle. The final recon- 
structed position (full circle) is a weighted average of those three points. The 
true source position is shown with a star. Open triangle shows the reconstructed 
source position for the top image if the head-tail discrimination failed. 



5.2. Relative light scale between both telescopes 

The mean amount of Cherenkov light on the ground pro- 
duced by a VHE gamma-ray shower depends mostly on its 
energy and the impact parameter (note however also the de- 
pendence on the relative position with respect to the magnetic 
North due to the geomagnet ic field effect, mostly pro nounced 



at lowest energies, see e.g. ICommichau et al.L 120081) . There 



fore, it is possible to investigate the relative light scale of both 
telescopes by selecting gamma-like events wi th a similar recon - 
structed impact parameter in both telescopes ( Hofmannl 2003 ). 
We apply rather strict cuts Hadronness < 0.2 and 8 2 < 0.01 to 
obtain a pure gamma-ray sample for this test. Hadronic back- 
ground events would spoil the resolution of this method due 
to the strong internal fluctuations and poor estimation of the 
impact parameter. We obtain that the absolute light scale for 
MAGIC II is ~ 10% larger than that of MAGIC I (see Fig. [8j. 
This difference is within the systematic error of the absolute 
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Figure 6: Theta 2 distribution of the signal (black points) in the analyzed sample 
of the Crab Nebula data and the background estimation (grey shaded area). 
Estimated energy 60GeV < E ea < lOOGeV (top panel) and E es , > 300GeV 
(bottom panel). The cut in theta 2 used for the calculation of the number of ON 
and OFF events is shown with the dashed, vertical lines. 



energy scale (see discussion in Section[6]l. 

5.3. Comparison of parameters 

We compare the rate of the reconstructed stereo background 
events above a given mean size (i.e. (Size\ + Size?)IT) for the 
data and the MC samples (see Fig. [9j. We normalize the rate o f 
the protons to the BESS measurement dYamamoto et al. , 2007 ). 
At the energy close to the MAGIC threshold, the rate of Helium 
nuclei in the cosmic ra ys (CR) is comp arable (~ 50 - 60%) 
to the rate of protons ( Horandel 2003I) . Helium nuclei can 
be roughly treated as 4 nearly independent protons of 4 times 
smaller energy. Thus the energy threshold for helium is higher 
than for protons and the rate of the MAGIC telescopes is dom- 
inated by the pro tons. Using proto n and Helium MCs, and the 
CR composition ( Horandel 2003 ) we estimate, that the effect 
of elements heavier than Helium on the rate of the background 
events registered by the MAGIC telescopes can be roughly ap- 
proximated by considering an additional 60% of Helium nuclei. 
After including all the components, the MC rate at low sizes 
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Figure 7: Distribution of the smallest charge of the brightest 3NN cluster for 
data and Monte Carlo (continuous and dotted lines respectively) for MAGIC I 
(top panel) and MAGIC II (bottom). 



agrees up to 20% with the data. Differences at higher sizes 
are caused mostly by the small systematic errors in the spectral 
slope reconstruction of both MAGIC and BESS and the limited 
maximum energy of the proton and the Helium MCs. 

The large data sample we use, and very good background 
reduction allows us to see a significant signal at medium ener- 
gies even with very loose cuts, which corresponds to a gamma- 
ray sample not biased by analysis cuts. Therefore, we are able 
to extract distributions of individual parameters of gamma-rays 
from the data and compare them with the MC predictions. The 
comparison of the Hillas parameters Width and Length of the 
gamma-rays excess computed from the data sample with those 
from MCs is shown in Fig. [10] for both telescopes separately. 
Those parameters mostly agree between the MC gamma-rays 
and the gamma-ray excess extracted from the data. 

5.4. Energy resolution 

We evaluate the performance of the energy reconstruction 
with the help of gamma-ray MCs. We fill histograms of 
(E rec - E, rue )/E, rue and fit them with a Gaussian distribution. 
The energy resolution can be calculated as the standard devi- 
ation of this distribution. In addition, the bias introduced by 
the method is estimated as the mean value of the distribution. 
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Figure 8: Relation between Sizei and Sizei for gamma-ray events obtained 
with a Crab Nebula sample. Only events with a difference in the reconstructed 
impacts of both telescopes smaller than 10 m are considered. Individual events 
with Hadronness < 0.2, 9~ < 0.01 and Impact parameters < 150 m are shown 
by black dots. The black, solid line shows the result of the fit, the dashed line 
corresponds to Size\ = Sizei- 



The dependence of the energy resolution and the bias on the 
true energy of the gamma-rays is shown in Fig.[TT] The energy 
resolution and bias weakly depends on the Hadronness and 6 2 
cuts, usually improving for stronger cuts. 

In the medium energy range (few hundred GeV) the energy 
resolution is as good as 16%. For higher energies it is slightly 
worse due to the large fraction of truncated images, and showers 
with high Impact parameters and worse statistics in the training 
sample. At low energies the energy resolution is worse, due to 
a lower photon number, higher relative noise, and worse esti- 
mation of the arrival direction, which spoils the precision of the 
Impact parameter reconstruction. In the medium energies the 
bias in the energy estimation is below a few per cent. At low 
energies (< 100 GeV) it is large due to the threshold effect. The 
bias is corrected in t he analysis of the s pectra by means of an 



unfolding procedure dAlbert et all 120071) . 



5.5. Angular resolution 

We define the angular resolution as the standard deviation 
of the 2-dimensional Gaussian fitted to the distribution of the 
reconstructed event direction of the gamma-ray excess. This 
corresponds to a radius containing 39% of the gamma-rays of 
a point like source. For completeness of the study we also cal- 
culate the 68% containment radius. In Fig.[T2]we show the Q 2 
distribution of the gamma-ray excess above two different en- 
ergy thresholds for data and MC simulations. The distributions 
for 2 < 0.025 can be reasonably well fitted with a single Gaus- 
sian. The distribution in a broader signal range is better de- 
scribed by a double Gaussian. In the lower energies MC and 
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Figure 9: Rate of stereo background events above a given mean size for data 
(the solid line) and MCs (dashed). The contribution of protons in the rate is 
shown with the dotted line, while helium and higher elements are shown with 
the dot-dashed line. 



data are consistent. For higher energies, due to small imperfec- 
tions not included in the MC simulations, the MC distribution 
is narrower resulting in < 20% higher number of events in the 
first bin (6 2 < 0.002). This does not influence the normal anal- 
ysis, since the usual cut in 8 2 is a factor of 5-10 larger, and the 
signal is recovered. 

The angular resolution obtained from MC and data are con- 
sistent (see Fig. [T3l , with just a small discrepancy mostly vis- 
ible at the highest energies (> 1 TeV). The stereo DISP RF 
method allows us to obtain a very good angular resolution of 
~ 0.07° at 300 GeV, and even better at higher energies. 

5.6. Spectrum of the Crab Nebula 

In order to minimize the systematic errors we apply 
Hadronness and 8 2 cuts with high gamma-ray efficiency for the 
reconstruction of spectrum of the Crab Nebula. The spectrum 
and the spectral energy distribution (SED) of the Crab Nebula 
obtained from the analyzed sample is presented in Fig. [14] The 
spectrum in the energy range 70 GeV - 1 1 TeV can be fitted 
with a curved power-law: 



dN_ 
~dE 



ME/300 GeV y +bl °s l0 (E/300 GeV > [cm -2 g -i Te y- 1 ] , ( 1 ) 



where the obtained parameters of the fit are: /o = (5.8+0.1 st at)x 
10- 10 , a = -2.32 + 0.02 stat , and b = -0.13 + 0.04 stat . 

The spectrum obtained with the MAGIC Stereo observations 
agrees within 20-30% with the previous measurements of the 
C rab Nebula. The curv ature seems to be smaller than reported 
in Albert et al. ( 2008bl) . however the effective spectral slope is 



still within the systematic and statistical error. No te moreover. 



that the fitting range is slightly different than in lAlbert et al 



d2008bl) . The spectrum in a broader energy range, together with 
its scientific implications will be discussed in another paper (in 
preparation), which exploits a dedicated analysis for the lowest 
and the highest energies. 
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Figure 10: Comparison of the Width (top panels) and the Length (bottom) parameters for MAGIC I (left panels) and MAGIC II (right). The solid line is computed 
with the MC gamma-rays, black points correspond to the gamma-rays excess, extracted with very loose cuts (Hadronness < 0.9, < 0.03) from the data sample. 
Only events with S ize > 200 phe are used. 



5.7. Integral sensitivity 

In order to allow for a fast reference and comparison with 
other experiments we calculate the sensitivity of the MAGIC 
telescopes according to a number of commonly used defini- 
tions. For a weak source and perfectly known background, the 
significance of a detection can be calculated with a simplified 
formula Access/ ^Nb g d, where N excess is the number of the ex- 
cess events, and Af Dg d is the estimation of the background. We 

as the flux of a source giving 



compute the sensitivity S 



xcess 

In the calculation of the sensitivity, S 



NexlyjNbkg 

5 after 50 h of effective observation time. 

, we are 



JNbkg. sys' 

applying also two additional conditions Access > 10, N excess > 
0.05A^b g d- The second condition protects against small system- 
atic discrepancies between the ON and the OFF, which may 
mimic a statistically significant signal if the residual back- 
ground rate is lar ge. Final l y, the sensitivity can also be cal- 
culated using the iLi & Mai (119831) . eq. 17 formula for signif- 
icance, which is the standard method in the VHE gamma-ray 
astronomy for the calculation of the significance. 

The integral sensitivity of the MAGIC telescopes working in 
stereo mode is shown in Fig. Q3] and the values both in Crab 
Units (C.U.) and in absolute units are summarized in Table [TJ 
Above a few hundred GeV the MAGIC Stereo sensitivity is a 



factor of 2 better than the one of the MAGIC I telescope. The 
improvement is even larger (by a factor of about 3) at lower en- 
ergies. Note that two identical telescopes working completely 
independently would have a sensitivity only a factor V2 bet- 
ter compared to a single one. The improvement by a factor of 
2-3 therefore comes from the efficient usage of the stereo pa- 
rameters in the analysis and the intrinsic reduction of the muon 
background in stereo systems. 

In Fig. [15] we also show the sensitivity obtained with sets of 
gamma-ray, proton, helium and electron MC simulations. We 
scale up the proton MCs by 20% to take into account the differ- 
ence in the rate between the data and the MCs (see Section [5~3l >. 
The flux of the electrons is much sm aller than the one of th e 
protons, and the spectrum is steeper (lAharonian et al. . l2009h . 



Nevertheless, due to their similarity to the gamma-ray showers, 
they can constitute a significant part of the background after 
gamma/hadron separation at medium and high energies (about 
15% at 200 GeV). The sensitivity obtained from the Crab data 
set is in agreement with the MC predictions. 

5.8. Differential sensitivity 

To evaluate the performance of the MAGIC telescopes for 
sources with a substantially different spectral shape compared 
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Table 1 : Integral sensitivity obtained with the Crab Nebula data sample above the energy threshold £thresh. ■ The sensitivity is calculated as Af cxccss / ^/Vbgd = 5 
after 50 h of effective time (S Nex i^j Nbkg ), or with additional conditions Access > 10. Nanus* > 0.05JVb e d ^ Nexl-^ Nbkg sys-*' ^ e sens itivity calculated to obtain 5cr 
significance after 50h of effective time according tojLi & Ma 1 1983|) with an assumption of 1 or 3 background regions is shown in Su&Ma.lOff an( ^ ^Li&Ma,30tf 
columns respectively. 
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Figure 11: Energy resolution (solid line) and the bias (dashed) obtained with 
the MC simulations of gamma-rays. Mild cuts are applied: Hadronness < 0.6, 
A 2 < 0.03. 



to the Crab Nebula, we compute also the differential sensitiv- 
ity, i.e. we investigate the sensitivity in narrow bins of en- 
ergy (5 bins per decade). In order to derive optimal cuts in 
Hadronness and 9 2 in each energy bin, we divide the sample 
into two roughly equal subsamples. In each energy bin we per- 
form a scan of cuts on the training subsample, and apply the 
best cuts to the second sample obtaining the sensitivity (see 
Fig.[T6]andTable|2]l. 

At medium energies, the sensitivity is very good (about 1 .5- 
2.5% C.U.). Even below 100 GeV, in the energy regime exclu- 
sive to the MAGIC telescopes among all the current IACTs, we 
obtain a good sensitivity of ~ 10% C.U. As in the case of the 
integral sensitivity, the differential sensitivity derived with the 
help of MC simulations is in agreement with the one obtained 
from the data sample. 
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Figure 12: Theta 2 distribution of the excess above 100 GeV (left panel) and 
above 300 GeV (right), for the Crab Nebula data sample (solid circles) and MC 
gamma simulations (open triangles). Fits with a 2D Gaussian distribution are 
shown with a solid (for data) or dashed (for MCs) line. A small 2D Gaussian 
mispointing with cr,,,,, = 0.014° was added to the simulations. 



5. 9. Off-axis performance 

Most of the observations of the MAGIC telescopes are per- 
formed in the wobble mode with an offset of 0.4°. However, 
serendipitous sources can occur in the FoV of MAGIC at an 
arbitrary angular o ffset from the viewing direction (see e.g. de- 
tection of IC 310, lAleksic et all 1201 Obi) . Therefore, we also 
study the angular resolution and the sensitivity of MAGIC at 
different offsets from the center of the FoV. Dedicated observa- 
tions of the Crab Nebula were performed at a wobble offset of 
<f=0.2°, 0.3°, 0.6°, 1°, and 1.4°. 

Since the presented here analysis is optimized for sources 
observed at the offset of 0.4° from the camera centre, the best 
angular resolution is obtained for this angle (see Fig. [171 . For 
larger offsets the angular resolution is slightly spoiled (~ 15% 
at £ = 1°), possibly due to higher influence of the optical aber- 
rations. MC simulations show a small (5% of the offset angle) 
radial bias in the estimated source position. Since normally at 
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Figure 13: Angular resolution as a function of the estimated energy obtained 
with the Crab Nebula data sample (points) and compared with the MC sim- 
ulations (lines) for the 2D Gaussian fit (solid line, full circles) and the 68% 
containment radius (dashed line, empty triangles). A small 2D Gaussian mis- 
pointing with o"„„, = 0.014° was added to the simulations. 
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Figure 14: Spectrum (top panel) and the SED (bottom) of the Crab Nebula 
obtained with the MAGIC Stereo syst em (black) compared with other experi- 
ments (gray curves): M AGIC I f solid. lAlbert et alll2008bl). HEGRA (dashed, 
lAharonian et alll2004l) . and H.E.S.S. (dotted. lAharonian et aljj2006l) . The ver- 
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Figure 15: Integral sensitivity of the MAGIC Stereo system, i.e. the flux of a 
source above a given energy for which Access / yf^bsd = 5 after 50 h of effec- 
tive observation time for MAGIC Stereo (black, solid - data, dashed - MCs) 
compared with the MAGIC I telescope (solid gray line). For comparison, frac- 
tions of the integral Crab Nebula spectrum are plotted with the thin, dashed, 
gray lines. 
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Figure 16: Differential sensitivity of the MAGIC Stereo system. We compute 
the flux of the source in a given energy range for which Access/ -JNbgi = 5, 
^excess > 10, Access > 0.05 A^gd after 50 h of effective time for the data (black 
points), and the MC (dashed line). The horizontal error bars represent the size 
of the bins in estimated energy (5 bins per decade). For comparison, fractions 
of the Crab Nebula spectrum are plotted with the thin dashed gray lines. 



tical error bars show statistical errors. 
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Table 2: Differential sensitivity obtained with the Crab Nebula data sam- 
ple. The sensitivity in each energy bin (E mm , £ max ) is calculated as 
^excess/ V^bgd = 5, Access > 10, Access > 0.05 Wbgd after 50 h. 
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Figure 17: Angular resolution (defined as the 68% containment radius) for 
gamma-rays with energy above 200 GeV for the observation of a source at a 
different offset from the center of the camera. Each point is obtained with a 
small, dedicated sample of Crab Nebula data. 



least two different pointing positions are used for each source, 
the impact of this bias on the data is reduced. 

The sensitivity above 290 GeV for different offsets is plot- 
ted in Fig. [18] The values of the sensitivity and the effective 
observation times of the data samples taken at different offsets 
are summarized in Table [3] The sensitivity of the MAGIC tele- 
scopes is quite good if the projected position of the source is 
within the inner part of the cameras. The sensitivity degrades 
significantly for sources located at a larger distance from the 
camera center (about a factor of 2 for 1° offset). 

6. Study of the systematic uncertainties 

Due to the complexity of the IACT technique, there are many 
factors which are only known with some uncertainties, thus 
contributing to systematic errors (see Table|4]i. We consider the 
systematic errors on the gamma-ray collection efficiency (i.e. 
on the absolute flux level), on the absolute light scale, and on 
the reconstructed spectral slope. 
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Figure 18: Integral sensitivity above 290 GeV for the observation of a source 
at a different offset from the center of the camera. Each point is obtained with 
a small, dedicated sample of Crab Nebula data. 





time[h] 




0.2 


0.46 


0.61 ± 0.15 


0.3 


0.23 


0.74 ± 0.25 


0.4 


8.7 


0.76 ± 0.03 


0.6 


0.89 


1.15 ± 0.14 


1 


1.55 


1.63 ±0.26 


1.4 


0.84 


2.8 ± 0.76 



Table 3: Integral sensitivity S 



N, 



above 290 GeV (calculated as 



Nex/-yNbkg 

:ess/ V^bed = ^ a ^ ter 50 h) obtained with the Crab Nebula data samples 



for different offsets f of the source from the center of the camera. 

6.1. Atmosphere 

The Earth's atmosphere has to be considered part of an IACT 
detector. Changing atmospheric conditions, small deviations 
of the density profile from the one assumed in simulations, as 
well as non-perfectly known atmospheric transmission due to 
Mie scattering introduce additional uncertainties in the energy 
scale and a small effect on the spectral slope reconstruction. 
We estimated the systematic error in the energy scale due to 
atmospheric parameters to be < 10%. 

6.2. Mirrors and night sky background 

An uncertainty in the amount of light focused by the mir- 
rors can be caused by non-perfect knowledge of the reflectivity 
of the mirrors, in particular short-term variations of the dust 
deposit. In addition to this, due to malfunctions of the AMC 
system, the total useful mirror area can vary from one night to 
another. Using the measurements of the reflected star images 
and the analysis of the muon images, we estimated that those 
effects produce a systematic error on the energy scale of about 
~ 8%. 

The level of the night sky background (NSB) can vary from 
one observation to another. In particular, galactic sources usu- 
ally have a higher NSB than extragalactic ones. Also, observa- 
tions in twilight and moonlight conditions exhibit higher NSB 
values. High NSB increases the fluctuations of the signal in a 
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Systematic effect 


Uncertainty 


F-Factor 


10% ES 


atmospheric transmission 


<10% ES 


mirror reflectivity 


8%ES 


PMT electron collection efficiency 


5%ES 


light collection in a Winston Cone 


5% ES 


PMT quantum efficiency 


4% ES 


signal extraction 


3%ES 


temperature dependence of gains 


2% ES 


charge flat-fielding 


2-8% ES FN 


analysis and MC discrepancies 


<10-15% FN 


background subtraction 


1-8% FN 


broken channels/pixels 


3% FN 


mispointing 


1-4% FN 


NSB 


1-4% FN 


trigger 


1%FN 


unfolding of energy spectra 


0.1 SL 


non-linearity of readout 


0.04 SL 



Table 4: Estimated values of the main sources of systematic uncertainties. 
They can affect the energy scale (ES), flux normalization (FN) and the spectral 
slope (SL). See text for the detailed explanation. 



pixel. Thus it can spoil the precision of the estimation of the 
Hillas parameters, and hence lower the acceptance for gamma- 
rays. We investigated this effect generating dedicated MCs with 
a 30% higher NSB than a typical galactic region of sky. Such 
an increase covers typical changes of the mean NSB observed 
in different regions of the sky in dark conditions, and a mild 
moonlight (moon phase < 20%). We found that the effect on 
the gamma-ray collection efficiency is up to 4% at low energies 
(~ 100 GeV) and negligible (< 1%) above 300 GeV. 

6.3. Cameras and readout 

An important systematic error in the absolute energy scale 
comes from uncertainties of the conversion coefficient of pho- 
tons to detectable photoelectrons. It contains uncertainty in the 
light collection in the Winston Cone (~ 4%), electron collec- 
tion efficiency of the first dynode (~ 5%), quantum efficiency 
of the PMT (~ 4%) and finally the F-Factor value of the PMTs 
(~ 10%). While it is difficult to disentangle and correct the in- 
dividual components of this en ergy scale uncertainty, by using 



an absolute muon calibration (Puhlhofe r et al.L 120031) . and in 



tertelescope cross-calibration we are able to obtain much more 
precise photon to phe conversion factors. 

In addition, a ~ 2% effect in the energy scale can be at- 
tributed to the temperature dependence of the gains. In order 
to make the response of the camera homogeneous, we perform 
the so-called flat-fielding. Flat-fielding equalizes the product 
of the quantum efficiency of the pixel for the wavelength of 
the light pulser with the gain of the PMT-FADC chain. The 
part of the signal coming to the trigger branch is thus only par- 
tially flatfielded, moreover temperature drifts can influence the 
flat-fielding. Those effects, and the fact that the flat-fielding of 
the signal chain is done only at one wavelength, can produce 
a ~ 6-8% systematic uncertainty in the energy scale and the 
event rates, mostly for the small events, and a smaller effect 
(< 2%) at higher energies. With over 1600 channels in total it 



is natural that a small number of them (< 10) at a given mo- 
ment is unusable, due to e.g. a malfunctioning PMT or readout 
electronics. Since we interpolate the signal in those channels in 
the analysis, the systematic effect on the energy scale is negli- 
gible. However the interpolation procedure can lead to a loss 
of about 3% of gamma-ray efficiency. Non-linearities in the 
analog signal chain and the small residual non-linearity of the 
DRS2 chip can produce a systematic uncertainty of about 0.04 
in the spectral index. 

6.4. Trigger 

The readout of MAGIC II introduces a dead time of about 
10% for typical trigger rate of events. The dead time is cor- 
rected in the analysis during the calculation of the effective ob- 
servation time by means of an exponential fit to the distribution 
of the time differences to the previous triggered event. The re- 
maining error is negligible (< 1%). Also, the inefficiencies of 
the trigger systems of the MAGIC telescopes are negligible. 

6.5. Signal extraction 

We investigate the systematic uncertainties induced by the 
signal extraction in MAGIC II by varying the size of the extrac- 
tion window. We find that the difference in the reconstructed 
number of photoelectrons is < 3%. It is similar to the corre- 
sponding value of MAGIC I. 

6.6. Mispointing 

Mispointing of one or both telescopes can influence the anal- 
ysis. Not only the 8 2 distribution becomes broader, but the rel- 
ative pointing differences between MAGIC I and MAGIC II 
spoil the reconstruction of the stereo parameters. The typical 
mispointing of the individual MAGIC te lescopes is < 0.02° 



dBretz et all l2009t lAleksic et all l2010ab . The final, post- 
analysis mispointing is slightly higher < 0.025° since it in- 
cludes reconstruction biases. We performed dedicated MC sim- 
ulations and conclude that the systematic error on the gamma- 
ray efficiency due to mispointing is < 4% 

6. 7. Background subtraction 

Due to dead pixels and stars in the field of view of the tele- 
scopes, the distribution of the events on the camera can become 
partially inhomogeneous (i.e. it loses its rotational symme- 
try). On top of this, the stereo trigger produced a natural in- 
homogeneity. As a first approximation, the stereo FoV can be 
treated as a non-circular overlapping region of two cones (size 
of each determined by the size of the camera trigger regions) 
originating from both telescopes. This results in an enhanced 
amount of events in the direction determined by the position of 
the second telescope and can lead to a 10-15% variation in the 
background estimation. Such an error would significantly influ- 
ence the reconstructed flux and spectral index of a weak source. 
However, since the MAGIC data are usually taken in wobble 
mode, alternating the source and anti-source position, the sys- 
tematic uncertainty of the background estimation is strongly re- 
duced down to < 2%. This procedure naturally does not work 
for serendipitous discoveries. Nevertheless, even in this case 



13 



a comparable decrease of the systematic error is obtained via 
renormalization of the source and anti-source (f plots in the 
off-source region or taking the background estimate from the 
geometrically equivalent position in the FoV of the other wob- 
ble position. The effect of inhomogeneity is mostly pronounced 
at the lowest energies, where usually the excess is smaller than 
the background. Thus, as an example, a weak source with an 
excess being just 25% of the background, the induced error of 
the flux estimation can be up to ~ 8%. 

6.8. Analysis 

Because of small remaining discrepancies between data and 
MCs the reconstructed spectra depend slightly (< 10% at 
medium energies, and < 15% at low energies) on the value 
of cuts used for extracting the gamma-ray signal. 

Also, the unfolding procedure, which is needed to correct for 
the finite energy resolution and the energy bias at the lowest en- 
ergies, introduces additional systematic uncertainties. We can 
estimate them by co mparing resul t s obta ined with different reg- 
ularization methods (I Albert et all 120071) . We obtain that for a 
typical power-law spectrum the spread of photon indices is 0.1. 

6.9. Total systematic uncertainty 

As some of the systematic errors depend strongly on the en- 
ergy of the gamma-ray showers, we summarize separately the 
systematic errors for low (< 100 GeV) and for medium en- 
ergies (> 300 GeV). As the individual errors are mostly in- 
dependent f r om eac h other we add them in quadrature (as in 



Albert et al. . 2008bh . The absolute muon calibration and in- 



tertelescope cross-calibration allows us to remove most of the 
errors connected with the photon to phe conversion. Thus, the 
energy scale of the MAGIC telescopes is determined with a pre- 
cision of about 17% at low energies and 15% at medium ener- 
gies. The systematic error of the slope of the energy spectra 
is estimated to be 0.15. At medium energies the error in the 
flux normalization (without the energy scale uncertainty) is es- 
timated to be 11%. At low energies the systematic errors are in 
general larger and the flux normalization is known with a preci- 
sion of about 19%. Additionally, we find that the Crab Nebula 
spectrum reconstructed by the MAGIC telescopes is consistent 
with other experiments within 20-30%. 

6.70. Run to run systematic uncertainty 

The total systematic error estimated in the previous section 
can be used when the data from the MAGIC telescopes are 
compared with the data of other instruments. However, large 
part of the studied here sources of the systematic errors, will 
result in a (nearly) constant offset which will affect all the data 
in the same way. Therefore, the remaining, relative systematic 
error which may change from one data run to another will be 
certainly smaller. In order to investigate this effect we divided 
our data into 40 min sub-samples. Such time binning allows 
us to compute the integrated flux (above 300 GeV) of the Crab 
Nebula with a precision of about 9% (see FigfT9l>. 

The standard deviation of the resulting distribution of fluxes 
computed from individual data runs is equal to 12%. While the 
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Figure 19: Integrated flux above 300 GeV of the Crab Nebula, as the function 
of the data sample number. The solid horizontal line shows the mean flux in the 
whole period. The dashed lines show the mean fluxes in November 2009 (left) 
and January 2010 (right). 



fluxes both in November 2009 and January 2011 can be well 
fitted with a hyphotesis of a constant flux (y 2 /«do/ = 8.8/8 
and 10.1/9 respectively), there is a clear shift of about 17% 
between both months. We conclude that the relative systematic 
error on short time scales (within a few days) with no drastically 
changing observational conditions is below 9%. On the longer 
time scales (months) the relative systematic may result in vari- 
ations of the estimated flux of the order of ~ 17%/ V2 = 12%, 
which is similar to the v alue obtained by the H.E.S.S telescopes 
( Aharonian et al. . 2006h . Note, that since the Crab Nebula may 
show intrinsic variability, the values derived in this section 
should be treated as an upper limit on the run to run system- 
atic uncertainty of the MAGIC telescopes. 



7. Conclusions 

We evaluated the performance of the MAGIC telescopes us- 
ing both MC simulations and a sample of 9 h of Crab Nebula 
data. We obtained that the energy threshold of the MAGIC 
telescopes (defined as the peak of the energy distribution for 
a source with a -2.6 spectral index) is 50 GeV. The upgrade 
of the MAGIC project with the second telescope led to a very 
good integral sensitivity of (0.76 + 0.03)% of the Crab Nebula 
flux in 50 h of the effective time in the medium energy range 
(> 290 GeV). An even bigger improvement (a factor of ~ 3 with 
respect to MAGIC I single-telescope observations) is obtained 
at lower energies, thus allowing to reduce the needed observa- 
tion time by a factor of 9. This makes the MAGIC telescopes 
an excellent instrument for observations of gamma rays with 
energies around 100 GeV. Thanks to the new Stereo DISP RF 
method, the angular resolution has improved as well (~ 0.07° at 
300 GeV). The energy resolution is as good as 16% at medium 
energies. 

We investigated different sources of systematic uncertainties 
and found that for a strong source like the Crab Nebula they 
dominate over the statistical errors. The spectrum of the Crab 
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Nebula obtained with the MAGIC Stereo system is consistent 
with the spectra of other VHE Cherenkov telescopes. 
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